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We investigate in details how the Virasoro algebra appears in the scaling limit of 
the simplest lattice models of XXZ or RSOS type. Our approach is straightforward but 
' to our knowledge had never been tried so far. We simply formulate a conjecture for 
the lattice stress-energy tensor motivated by the exact derivation of lattice global Ward 
identities. We then check that the proper algebraic relations are obeyed in the scaling 
limit. The latter is under reasonable control thanks to the Bethe-ansatz solution. The 
results, which are mostly numerical for technical reasons, are remarkably precise. They 
are also corroborated by exact pieces of information from various sources, in particular 
Temperley-Lieb algebra representation theory. Most features of the Virasoro algebra (like 
central term, null vectors, metric properties...) can thus be observed using the lattice 
models. This seems of general interest for lattice field theory, and also more specifically for 
finding relations between conformal invariance and lattice integrability, since basis for the 
irreducible representations of the Virasoro algebra should now follow (at least in principle) 
from Bethe-ansatz computations. 
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1. Introduction 



In the last years there has been much progress in the study of two dimensional inte- 
grable lattice models and integrable field theories. Deep analogies have been discovered 
between the two subjects. 

It is certainly expected that some relations exist between integrable lattice models and 
integrable field theories because of their many common mathematical ingredients. Such 
relations are often straightforward; for instance R matrices appear both in the construc- 
tion of Boltzmann weights []T| and the study of braiding properties of conformal fields 
simply because of the ubiquity of the Yang-Baxter equation. There are however deeper 
relations that are not yet understood. For instance in § it was shown that local height 
probabilities (or order parameters) in restricted solid on solid (RSOS) integrable lattice 
models can generally be expressed in terms of branching functions of conformal field theo- 
ries. A consequence of this observation is that the logarithm of the corner transfer matrix 
^ provides a "lattice Virasoro generator Lq" . More precisely observe first that the natural 
configuration space of RSOS models is a space of paths on restricted weight lattices 0. 
For any representation of the Virasoro algebra (we restrict to this algebra for simplicity) 
there exists a choice of boundary conditions for the RSOS paths such that, restricted to 
these paths, the trace of the exponential of the corner transfer matrix coincides with the 
Virasoro character. This may be some sort of coincidence due to the common presence of 
elliptic functions in the parametrization of off-critical integrable Boltzmann weights and 
in the branching functions. However the final formula for Lq and the path structure of 
the potential Virasoro representations are so beautiful that one hopes there exists in fact 
a complete action of the Virasoro algebra on the space of paths. This would be physically 
intriguing because the Virasoro algebra is traditionnally associated with conformal invari- 
ance, a property of critical continuum theories, while we deal here with discrete, off-critical 
models. This would also be mathematically useful as the RSOS paths probably would pro- 
vide very simple basis for the irreducible Virasoro modules and maybe play a role similar 
to the crystal basis The action of the Virasoro algebra on RSOS paths has only been 
exhibited in the Ising model using the underlying fermions [§ so far, although there 
has been much progress recently in understanding the spectrum of corner transfer matrices 
||10|| ||1 1|| . The possible action of a current algebra on paths has been studied in |T^ . 

Close relations between integrable lattice models and conformal field theories have also 
been observed is the spectrum of row to row transfer matrices in critical models. In [0| 



(see also [^) it was found for instance that degeneracies in the free field representations 
of the Virasoro algebra could be observed as well on finite lattice systems of XXZ type. 
Qualitatively, the lattice quantum group plays a role similar to the screening operators 
algebra and its commutant, the Temperley-Lieb algebra, a role similar to the Virasoro 
algebra. So far, only Virasoro characters have been reproduced, by taking a 'scaling limit" 
(to be discussed in details later) and using for "lattice Lq" the logarithm of the row to row 
transfer matrix. The situation is thus rather similar to the one we discussed previously for 
off-critical models and corner transfer matrices, and the most intriguing question is again 
to find some expressions for the other Virasoro generators. A first difference is that the 
results obtained from the row to row transfer matrix are not as nice as the ones ones from 
the corner transfer matrix. Another difference is that we deal now with a critical system, 
where appearance of the Virasoro algebra in the scaling limit is expected. The search for 
the Virasoro algebra in the lattice model ITB] gains thus important physical motivations. 
Some natural questions to ask are: is there some nice lattice deformation of the Virasoro 
algebra? is there a reformulation of lattice integrability as some sort of lattice conformal 
invariance? are there lattice Ward identities? 

In this work we investigate in details how the Virasoro algebra appears in the scaling 
limit of the simplest lattice models of XXZ or RSOS type, and give a partial answer to the 
third question. For a different approach related to the same problem see [jl^. The paper 
is organized as follows. In section 2 we discuss the concept of lattice stress energy tensor. 
By generalizing the ideas of Kadanoff and Ceva we study the behaviour of the square 
lattice Q-state Potts model under stretching and straining and we write global lattice 
Ward identities. This allows us to make a reasonable guess for the local stress energy 
tensor. We reformulate this algebraically using the Temperley-Lieb algebra formalism. 

The strategy is then to simply compute commutators of our lattice expressions and 
study their behaviour for large systems. 

Section 3 contains first the detailed example of the Ising model where as usual things 
work out nicely. We then discuss in some generality what kind of convergence can be 
expected for the lattice quantities, and we define precisely the scaling limit that should 
be taken. Sections 4 and 5 contain a detailed study of the XXZ chain. In that case the 
lattice model behaves like the Coulomb gas with a charge at infinity and provides non- 
unitary representations of the Virasoro algebra. We check numerically that our lattice 
expressions disentangle the Left and Right sectors and reproduce the necessary commu- 
tation relations, in particular the central term and the null vectors structure. Section 6 



contains a similar study for the periodic RSOS models. In that case the Temperley-Lieb 
algebra representation is unitary and so is the Virasoro representation obtained. Section 7 
contains finally a discussion of the case of fixed boundary conditions for the RSOS models. 
Technically this is the most difficult case because convergence of numerical data is not as 
good as for periodic boundary conditions. It is however far better from a conceptual point 
of view since in that case there is a single Virasoro algebra and by appropriate choice of 
boundary conditions a single irreducible representation can be selected. Section 8 contains 
final comments, in particular about the possibility of making the scaling limit simpler by 
considering time dependent correlation functions. 



2. Lattice stress energy tensor 

We now discuss the concept of lattice stress energy tensor for critical integrable lattice 
models based on the Temperley Lieb algebra. Our arguments are a simple extension of 



the pioneering work of Kadanoff and Ceva [18 



2.1. Q state Potts model: generalities 

We consider the Q-state Potts model on a square lattice whose vertices are labelled 
by a pair of integers (j, k). To each vertex is thus associated a spin variable ajk that takes 
values 1, . . . , Q. The exponential statistical weight has the form 

A = Y, K^ij + 1/2, k)6{ajk, + Ky{j, k + l/2)6{ajk, a^, fc+i), (2.1) 

jk 

where the couplings are in general edge dependent and the natural labelling of edges is 
illustrated in figure la. The usual partition function is 

^(W) = E^^' (2.2) 

CTjk 

where the notation indicates dependence on all the couplings and the sum is taken over all 
possible values for each spin. We introduce a more properly normalized partition function 

as 

n{K})= n (^)'''^' (2.3) 
edges 

where the notation means that a product over all edges {j + 1/2, k) and (j, k + 1/2) is 
taken and for each, the proper coupling constant put in denominator. By standard high 



temperature expansion techniques the Q-state Potts model can be defined for any real Q 
i ■ 

An important property is the duality symmetry. Define the dual lattice as the square 
lattice whose vertex (j, k) stands in the middle of the face of the original lattice with south 
west corner (j, k) (see figure lb). Define the mapping K K*{K) by 

^ (24) 

Define accordingly couplings on the dual lattice 

K*ij + 1/2, k) = K'iKyU + l,k + 1/2)), (2.5) 

and 

K;{j, k + 1/2) = K*{K,{j + 1/2, k + 1)). (2.6) 
The dual of the original Potts model is the model with couplings ( |2.5| ) , ( pj.6D defined on 



the dual lattice. As follows from the standard analysis |^ the partition functions for the 
original model and its dual are simply related: for large systems where boundary effects 
can be neglected one has 

Y{{K})=Y{{K*iK)}). (2.7) 

We shall call homogeneous the case where couplings depend on x,y but not on the 
position. It is convenient to introduce the variables 

_ e^- - 1 _e^y -1 



One has then 

Y{{K}) ^ y(e., ty) = Y{ey, e.) = y (1 1) , (2.9) 

where the last two equalities hold for a big system and follow respectively from x ^ y 
symmetry and duality. 

In the homogeneous case the model is critical when Q 

= K*{Ky), Ky = K*{K^). (2.10) 

We are interested only in second order phase transitions which occur for < Q < 4. It is 
convenient to parametrize then 

^ = 2cos7,7G [0,7r/2], (2.11) 

and 



^» = - - -zP^, (2.12) 



1 _ sinw 
sin(7 - u) 

where u is called spectral parameter. We shall restrict to the ferromagnetic regime 
Kr,,Ky> for which u G [0, 7]. 



2.2. Lattice T^x 

When the model is at a second order phase transition point, all correlation functions 
at large distance depend on a single variable with elliptic like symmetry 

7^2 = +52^2_ ^2.13) 
The parameter S can be determined in various ways. To start, we simply borrow the result 



of PI 

= cotan ( —) . (2.14) 



.27; 

As a result, a derivative with respect to S is equivalent to straining the system since for 
functions that depend on ( |2.13D one has 

d d d 

Refer now to the relative roles of temperature and energy for analogy: the integral over 
the system of the energy density is an "operator" which acts (when inserted in correlation 
functions) like a derivative with respect to temperature. Look here for an operator which 
act s as 

(2.16) 

and similarly for any local quantity other than the spin cr. In this section < . > indicates 
thermal average. Relation ( |2.16|) is expected to hold in the limit where j, /c >> 1 so these 
variables can be treated as continuous. From ( p.l5|) we can perform the right hand side 
operation in ( |2.16|) simply by taking instead derivatives with respect to S. From ( |2.14| ) 
this is equivalent to taking derivatives with respect to u, and from ( |2.8D and ( |2.12| ) this is 
also equivalent to taking derivatives with respect to K^, Ky. But such derivatives can be 
taken by acting with an operator that depends on the spins only. Explicitely things are as 
follows: one finds first 



9 47 1 sin 7 



S— = 



e^- dKx e^y dKy 



(2.17) 



dS 71(52+5 2 ginMsin(7 — m) 

Now to take derivatives with respect to K^, Ky one just has to insert 5 functions of the 
appropriate spin interaction: for instance 
d 

dK {j + 1/2 k) ^ >= [< S{(7jk,crj+i^k)- < 5{ajk,(7j+i^k) >] f^iifci^j2fc2 

(2.18) 



and similarly for any local quantity. Thus introducing 

47 1 sin 7 



h{j + l/2,k + l/2) = 



TT iS^ + 5 2 sin u sin(7 — u) 



0+1/2, fc) _ I ^Ky(j,k+i/2) _ I 

(defined up to an additive constant) one has 



(2.19) 



Ti = 5^ti(j + l/2,A: + l/2). (2.20) 

jk 

Obviously Ti changes sign in a; ^ y. On the other hand it is invariant in duality transfor- 
mations since 

dK* e-^* - 1 



dK e^* - 1 ■ 

We now want to show that equation ( |2.19|) , which can easily be guessed using the 
expected behaviour of 71 in x ^ y and duality transformations, leads to the result ( |2.14|) . 
Observe to start that the homogeneous lattice model with coupling constants Kr^, Ky can 
be considered at large distance as discretization of an isotropic continuum by rectangles of 
unit area, whose shape can be characterized by an anisotropy angle 9 as shown on figure 
2. Clearly one has 



= tan y-j . (2.21) 

The problem is to determine the dependence 9{Kjc,Ky) on the self dual line ( |2.1C1|) . To 
do so, recall from the definition of the stress energy tensor in a continuum theory (we use 
the normalization which is standard in conformal field theory [^) that 



1 f / d d \ 

^^^^ ^2^ J V dx^^ ^ d^^'^ / ^^'^^^ 

for a deformation x ^ x + 5x and y y + 5y. Considering the particular deformation 
^ ^ ^ + 5^ one finds 

d l cos^/2 ^^ d 1 sing/2 

a^'^" = 2^1^'^^' d-y^y = -2^2^^- ^2.23) 

Hence, using also Txx "l~ — 0, 

5lnZ If , , . . 

-TT- = TT— ^ <Tx.> dxdy. 2.24 
dtl An sm 9 J 



For the moment simply observe that the variation of In Z is proportional to the integral of 
< T^^ > which is itself proportional to Ti (this follows from the definition (eg ( p.l6|) ; see 



also the subsequent discussion of lattive versus conformal Ward identities). On the other 
hand, since 9 can be expressed in terms of K^;, Ky on the self dual line, and assuming this 
expression is invertible, one has 

where 

dK^ dKy 

Now comparing ( p.25|) , ( p.24|) and ( p.l9| ) we see that the dependence 9{Kx,Ky) must be 
such that 

G — 1 

It is easy to show that this implies, from parametrization ( |2.12| ) that 9 is proportional to 
7 — M. To determine the proportionality coefficient simply observe from ( |2.12| ) that the 
singular values 9 — 0,9 — ir must correspond to u = 0, u = 7. Hence 

e=^lh^. (2.28) 

7 

2. 3. Comparison of lattice and conformal Ward identities 



We now compare the lattice Ward identity ( p.l6 ) with local Ward identities in con- 
formal field theory. Recall the definitions 

T{^z) = — {Txx ^yy '^^Txy) , (2.29) 

and 

f{z) = ^ {Txx - Tyy + 2iTxy) , (2.30) 

and the local Ward identity 

{Txx{z, Z)(t)i{zi, Zl) . . . , (t)N{zN, Zn)) = 

^ (I^ + (i^ + ai- + ai- j • • ■ ' ■ 

(2.31) 



Now integrate ( |2.31j ) over the whole complex plane. Using 



Z - Zi 



one can -formally- establish that 



dxdy 



Z - Zi 



while by derivation 



Now use 



dxdy 
{z - ZiY 



-TlZi, 



0. 



(2.32) 



(2.33) 



_ d ^ d d d 
dz dz dx ^ dy' 



to get, integrating (|2.31| ) 



i=N 



T^x<Pi ■ ■ ■ (pNdxdy"^ ^ ^ ^ (^~^*^ 



(2.34) 



We could also have obtained this relation directly from the classical definition of the stress- 
energy tensor. This looks like equation (|2.16|) with the correspondence 



T^-<T^>^ / ^^^^^^dxdy. 



(2.35) 



The only condition for the lattice Ward identity to hold is j, /c >> 1. 

So far, only 7i, a global quantity, has been identified. We do not know any convincing 
way to establish a priori that the density ti ( |2.19|) is indeed what will correspond to the 
local quantity T^^. We shall assume it is true for the moment, that is 



h{j + 1/2, k + 1/2)- < h >^ — (T,, - Tyy] 



(2.36) 



Most of the paper will be devoted to checks of that hypothesis. The precise meaning of 



the correspondence indicated by i-^ is not as simple as for the global result (|2.35|) . This 
leads to the proper definition of the scaling limit to be discussed in the next section. The 
same remarks applies in all the rest of this section for expressions involving the notation 



2.4- Temperley-Lieb algebra and operator formulation 

To check that ti corresponds to Tr^j; a first possibiHty would be to investigate directly 
local Ward identities for the lattice model. This however requires the knowledge of many 
multipoint lattice correlation functions and is therefore very difficult, analytically or nu- 
merically. It is much easier to study the algebraic relations satisfied by the modes of ti 
and compare them with the Virasoro commutation relations. To do so however we need 
to pass to a hamiltonian description of the lattice model. 

Consider the row to row transfer matrix for the Potts model, with time in the x 
direction. The transfer matrix is a true operator acting on configurations of vertical spins, 
which we call (discarding now the horizontal index) aj, j = 1, . . . , L where L is the size 
of the system. We choose periodic boundary conditions in space direction. Introduce the 
operators 

X2j = 1 + eye2j, X2j-i = \fQ (e^ + e2j-i) , (2.37) 
where are defined by their matrix elements 

(e2j)^^^, = v^Jj5(afc,4)5(aj-,(jj+i), (2.38) 

k 

(with CTj^-i-i = o\) and 

(e2,-i),,, = ^n'^(^^'^^)- (2.39) 

They satisfy the relations 

= ejej±iej = e^, [e^-, Cfc] = 0, \j -k\>2, (2.40) 

where e2L+i = ei. The transfer matrix itself reads 

j=L j=L 

^= n^2, nx2,-i. (2.41) 

3=1 i=i 

In this formalism, ti takes a rather complicated form in general. It simplifies however in 
the hamiltonian limit where Ky and Kr^ oo. This corresponds to w ^ 0. In this 
limit ey ~ ^y^- K is then convenient to proceed as follows. Call F the prefactor in ( |2.19| ) 
and (^) . Then 



(ti(ii + 1/2, h + l/2)ti{j2 + 1/2, k2 + 1/2)) = 



e^- -1 d 



-1 d 
e^y die, 



X 



-1 d 



-1 d 

e^y dK. 



In Z, 



(2.42) 



y J 



where the labels of the couplings K are implicit. For a general heterogeneous model the 
partition function on a rectangle of size L, T with periodic boundary conditions in time 
direction as well reads 

Z = tr (f X . . . X f) , 

where the product is over T terms and all the couplings dependence is implicit. Compute 
( p.42| ) and identify the result with the correlation function of two operators introduced at 
the appropriate places. Using the relations (|2.40| ) one finds 



ti ~ — rey(e2j + e.2j-i) + constant. 



(2.43) 



To compare with the operator Tn^j^ a final rescaling has to be performed since the units of 
length in x and y direction are different: 



cotan ^0 [ii- < ii >) ^ ^ {Txx - Tyy 



(2.44) 



From the expression of F and ( ^.361 ) one finds the central result 



27 1 ^ 

: {e2j + e2j-i - 2eoo) ^ t^T'xi;, 

7rsm7 ^ ^ 2tt 



(2.45) 



where e^o is the mean value of the e operators in the ground state for L ^ 00. This value 
is non universal. Let us emphasize that in the conventions used so far, the system has 
length L. We will change our conventions to 2L later. 



2.5. Lattice 



xy 

So far we have studied only one component of the stress energy tensor. The other one 
should satisfy 

f d d d d \ 

< {T2- < T2 >) aj,k,CTj,k, >= U'l^ + h-g-+j2-gj^ + k2-^ 1 < cTj.k.crj.k, >, 

(2.46) 

and be even under duality and x ^ y transformations. A plausible expression that satisifes 
these constraints is 



t2ij + l/2,k) oc 



^K^{j + l/2,k) _ I ^Kyij + l,k+l/2) _ I 



d 



d 



gK^{j+i/2,k) gX,o+i,fc+i/2) dK^{j + 1/2, k)dKy{j + l,k + 1/2) 



^K^{j + l/2,k) _ ^ gKj,(j + l,fc-l/2) _ ^ 



d 



d 



gK4i+i/2,fc) gK,o+i,fc-i/2) dK^{j + 1/2, k)dKy{j + l,k- 1/2) \' ^^■^'^^ 
where it is necessary to introduce two derivatives to ensure x ^ y invariance. Another 
plausible expression would be 

"e^x(i+i/2,fc) _ I ^Ky{j,k+i/2) _ Y 



t2{j + l/2,k) (X 



d 



d 



gK.(i+i/2,fc) gK^(i,fc+i/2) dK^{j + 1/2, k)dKy{j,k+ 1/2) 

gi<r,(j+l/2,fc) _ I ^Ky{j,k-l/2) _l Q Q 



gK.0+i/2,fc) gK,o-,fc-i/2) Qx^{j + 1/2, A;) A; - 1/2) 

. Corresponding to ( |2.35| ) we should have 



To- < T2 >v 



T +T 



27r 



dxdy, 



and, assuming the correspondence extends locally 



t2{j + l/2,k)~<t2 >^ ^{T^y+Ty^). 



(2.48) 



(2.49) 



(2.50) 



An formula similar to (p.43| ) can be derived for these two expressions which from 
( ^.38| ) and ( ^.39] ) involves commutators [e2j, 62^-1] and [e2j, 62^+1]. The result is 

2 



7 



Ztt 



(2.51) 



([e2j-i,e2j] + [e2j,e2j+i]) 

^7rsm7/ 

(by reflection symmetry there is no mean value to subtract). 

2.6. Conserved quantities 

For the homogeneous model on the self dual line, the row to row transfer matrices 
f{u) do not form a commuting family . To obtain a commuting family it is necessary, 
like in the Ising case , to consider diagonal to diagonal transfer matrices (see flgure 
3) 



Td - Yl^2jX2j + l. 



(2.52) 



The matrix f^iu) can then be considered as a generating function for an inflnite family of 
commuting hamiltonians. The flrst few follow from 



In [Q-'^/h^TDiu] 



u 



2L 



2L 



sm 7 ^ 2 sin 7 ^ 



2L 



U 



3 sin 7 
Defining generally 



oo / \ 1 

^ \ _L/,(n) 



n=l 



sm 7 / n 



(2.53) 



^2.54) 



it is possible to get all the h^"^^ by using a ladder operator. Introduce 



•J' 



(2.55) 



then the generic term of h^^^ is obtained by computing the generic term in [/i('^-i),L]. This 
follows from a simple generalization of the arguments in (the statement has indeed to 
be given for generic terms because there are boundary effects for a finite chain). 

One the other hand we have already explained how the spectral parameter is related 
to the anisotropy angle 9 (|2.14 ) in the continuum ( p.21| ) . As a result, consider the transfer 
matrix Vtd where V is the translation operator of half the diagonal of an elementary 
rectangle. For a conformal field theory with the same geometry (figure 4) we would have 
the evolution operator (at imaginary time) 



U = exp 



2 

2^ 



sin 9Tyy + cos OT^v ] dv 



(2.56) 



where we have rescaled by a factor 2 due to the diagonal geometry, u, v are axis rotated 
by 7r/4 with respect to x,y. Therefore we expect 



7^fi:»(M) ^^ exp sin6' ^Lo + Lo - + zcos6' (Lo - Lo) |. 

Expanding and comparing with ( |2.54|) we have 



^2.57) 



2n-l) 



and 



2n) 



TT sm7 

7 



7rsm7 



2n-l 



TT 



IV- Lo + Lo 



12 



2ti 



-1 



TT 



7 y L 
These results agree in the case n = 1 with ( |2.45|) and (|2.51| ) 



(2.58) 



(2.59) 



2.1. Remarks 

Let us just notice that 



£(j + l/2,A;)cx 



^K^{j+l/2,k) _ Y ^Ky{j,k+l/2) _ I 



(2.60) 



is invariant in x ^ y and odd in duality and is a reasonable candidate for the energy 
operator in the Q state Potts model, e being odd in duality implies that all its n point 
functions with n odd vanish, a result in agreement with the known result that correlation 
functions of e in the lattice model coincide at large distance with those of the operator 
]T3| in the corresponding conformal field theory (where (prs denotes the primary field 



-'12 



of conformal weight hrs)- 

2. 8. Other representations of the Temperley-Lieh algebra 

Besides the Potts model, several other statistical models are obtained by choosing 
different representations of the Temperley-Lieb algebra. Such models include the 6 vertex 



model [H , the loop models |21[] , the ADE restricted solid on solid models [^. The oper- 
ator expressions like ( |2.43| ) can still be defined in these other models, and we believe that 
they correspond to their stress energy tensors. In fact choosing different representations 
of the Temperley Lieb algebra will simply correspond to choosing different representations 
of the Virasoro algebra. 



3. The scaling limit: generalities 

In this section we discuss in what sense ti can be considered a lattice approximation 
for Tn^^ and define a scaling limit into which the Virasoro commutation relations will be 
observed, hence giving a meaning to the notation h-> of the previous section. 

3.1. The Ising model 

The representation of the Temperley Lieb algebra corresponding to the Ising model is 
obtained by the quotient ||2^ [Q] 



ejej+i + ej+iCj - v2(ej + Cj+i) + 1 = 0, (3.1) 
which can be parametrized by 

e, = -^+zv^r,r,+i, (3.2) 



where 



{r.,r,} = 5. 



(3.3) 



The periodic RSOS model corresponds to antiperiodic fermions. It is easy to see then that 
aU the in ( 2.53| ) are hnear combinations of terms of the form T jV jj^2p+i with p < n 

while the /i*^^") are combinations of terms of the form TjTj-^.2p, with p < n. It is therefore 
enough to study the hamiltonians 



H 



2L 



One has in particular 



and 



2 2V2 



8i 

We shall more generally consider the quantities 

. 2L 

where / is a C°° function of period one and we have used the short hand notation 

p- 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



/ 



J + 



f 



2L 



(3.8) 



The function / has Fourier representation 



f(x) = J2fie 



2i7rZa 



(3.9) 



while we set 



with 



and thus 



2L f 



(3.10) 

(3.11) 
(3.12) 



By elementary computation one finds 



=H — <u, H u> 



J2 i-ir sin {^{k' - k)) : vCfcVCfc, : 6{k + k' + l, 2qL) 



\2L 

k,k' ,q 



(3.13) 

with 

< u,H^"^u >= -^5(p odd) . \ (3.14) 
where the normal order is defined by 

: : = -^j-^k: k > j 

(3.15) 

= -^k'^j, k<j 



and the ground state |u > is annihilated by with k > 0. In this section the metric is 
such that |u > has norm one and = ^E'-fc (coinciding with the natural metric for Ising 
spins). One has 

<„,r,r.u>=-i-^y-^. (3.16) 

In this case, what we want to call the scaling limit is easy to understand. First one 
restricts to states |v > such that 

^'fc|v >= 0, k = ko + ^,...,L-ko-^, (3.17) 

where ko is kept fixed as L ^ oo. Suppose now we choose the function / to have a single 
Fourier component : /; = 6i^n- Then one finds, by keeping also n fixed as L ^ oo 

jj(p)^^2™) ^I_J2(^k + ^) : a-kak+n : -(-1)^ " i) ^ ^' (3-18) 

k 

where we have set 

(3.19) 

Due to the restriction ( |3.17| ) there is no ambiguity in this definition for L large enough. 
After taking the limit L — > cxd one takes the limit k^ ^ oo to obtain indeed the well known 
representation of the Virasoro algebra in the Neveu Schwartz sector (see eg ||2^) that is 



Hf^'\e2™)^^(L„ + L_„) 
^e^™) I (L. - L_) 



the excitations near k = and k = L corresponding respectively to the right and left 
sectors. It is this double limit process that we call scaling limit in the following. 
Consider now the commutator 



-/ 



i=i 



S 77,2 

2 ~ YJ 



9 



s ni 



+ f 
-f 



S 7l2 

d n2 



9 J + 



ni 



j + 



d 77-2 



d ni 
d ni 



r r 



r r 

J- j-L j+n2-ni 



where we have set 

77,2 + ni = s 
n2 — ni = d 

Expanding the functions to first non trivial order we find 

2L 



(3.21) 
(3.22) 
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j + 



j + 



ni 



/ 



j + 



j + 



9 



9 



j + 



j + 



^2 ( - <U, Tj-r^+^^+^jW >) + 



(3.23) 



ni + r 



' . d' 




' . d 




9 





n2 



j+n2 — ni ■ 



< u,TjTj+^,.^,u >) + ... + <u, [Ri^'\f),Ri^'\9)]u > 



where dots stand for higher order terms in the expansion of the functions. For each of 
the two terms in the right hand side of the above sum we can apply the previous analysis. 
Replace each sum by its (exact) representation ( |3.13D and then consider the scaling limit 
for the right hand side as in ( |3.18| ). One gets then 



71 TC 



1 L 

2m TT 



-H(-)(/),-Hi-)(<7) 

TT TT 



(3.24) 



u > 



The first term on the right hand side appears only in the scaling limit: for finite L the 

- (n) 

algebra of lattice quantities (/) does not close. So far we have kept L finite in the 
second term. Dots stand for the effect of higher order terms in the expansion of the 



functions, and corrections to the scaling limit in ( |3.13|) , both being negligible. Averages 
of product of gamma operators are computed using ( |3.16| ). We look for all contributions 
that do not vanish as L ^ oo. It is easily seen that there are of two types. Keeping only 
the leading order in ( p.l6|) produces a term with third order derivatives of functions /, 



all previous order giving vanishing contributions by periodicity. Keeping the next order 
in (|3.16|) produces a term with first order derivatives of functions /, g. These two types of 



terms are of order L^. All subsequent contributions vanish as L — >^ oo. One finds therefore 




TX TT 



u > -- 

2L 1 „ /•2L 



(2z7r)3 12 



/ {f"'g-fg"')dv-:^^ [ {f g' - f'g)dv] 6{m + n2 odd) + 0{L-^) 
Jo 1^ Jo J 

(3.25) 



3.2. Scaling limit 

Let us summarize our observations in the Ising case. A first difficulty when one deals 
with a lattice model is to "follow states" for increasing system sizes to be able to define 
a limit process. This is easily done in the Ising model if one characterizes states by the 
fermionic modes which are occupied. In that case one can also, by identifying the lattice 
and continuum fermions, know which lattice state goes to which continuum state in the 

(n) 

large L limit. We then have lattice operators like h^'^^ (/) or equivalently the H (/) which 
have the following behaviour. If one first restricts their action to a finite number of low 
lying excited states characterized by some parameter ko ( |3.17|) and let L 00 the action 
of these operators on the selected states coincides with the action of Virasoro generators on 
the states of the continuum theory which have been identified as the limits (in the above 
sense) of these lattice states. One can then let ko go to infinity to recover the complete 
Virasoro action ( p.l8|) . This process is what we call scaling limit of lattice operators. Of 



course the algebra of lattice quantities does not close for finite L. However, in the Ising 
case, it closes once the scaling limit is taken for the right hand side too. Moreover the 
central term can be computed by taking the L —>■ 00 limit of the central term obtained in 
the lattice commutator. 

The definition of the scaling limit is easy to generalize in principle. This is done in 
practice using the Bethe-ansatz. The latter provides a way of organizing lattice states and 
identifying them through sets of integers which bear some resemblance with the fermionic 
modes of the Ising case. We will then exhibit lattice operators whose scaling limit will be 



the Virasoro generators. However one has to be careful in trying to close the algebra of lat- 
tice quantities, because the scaling limit of a commutator is not in general the commutator 
of the scaling limits. To close the algebra, we will need to restrict also the intermediate 
states in the computation of commutators. 

We illustrate this point in the next subsection, before embarking on the study of the 
scaling limit for lattice Virasoro generators in the next sections. 

3.3. Scaling limits and commutators 

Suppose that we compute the complete commutator of lattice quantities using the 
lattice stress energy tensor . How wrong would the result be? we show in this section that 
for the measure of the central charge, it would in fact be surprisingly acciirate for x not 
too large. 

Consider therefore the hamiltonian h^^^ (|2.58|) and set more generally 



2L 



Similarly introduce 

2 2L 



P(/)--(^j E/b + l/2][e„e,+,]. (3.27) 



Consider now 



m)M9)] = (-^) ' J2(f^M + l]-f[j + MMej. e,+i]. (3.28) 
\ TT sm 7 / ^—^ 



By expanding the combination of functions /, g one finds 



-H(/),-H(^7) 

TT TT 



Consider now 

3 2L 



^ ^P{f'g-fg')+0{L-'). (3.29) 



21% TT 



1 / \ . , 

[H(/), P(^7)] = -- ( ^ j - 3/2][e„ [e,_2, e,_i]] 

+g[j - 1/2] [e„ [e,_i, e,]] + g[j + 1/2] [e ]] + g[j + 3/2][e^,[e^+i,ej+2]]} 



(3.30) 



Let us investigate the vacuum expectation value of this quantity, denoting by |u > the 
ground state. The metric is such that |u > has norm one and is for instance the nat- 
ural metric in the vertex model representation (see later). One finds, using translation 
invariance and the Temperley Lieb defining relations 

<u, [^(/),P(^7)]u>= 

3 



1 



7 



2L 



<n,[ek,[ek+i,ek+2]]vi>y2f[j]{g[j + 3/2] - g[j -3/2]) 

2L 

+ < u, (^v^(efcefc+i -fefc+iCfc) -2efc) ,u> ^ + 1/2] - - 1/2]) 



Consider now the case 



(3.32) 



so 



H(/)^^(Lp + L_p) 



(3.33) 



If the scaling limit of the commutator was the commutator of the scaling limit, the right 
hand side of ( |3.31j ) would be, at dominant order 



TT 



Jp+q 



2p < u, (Lo + Lo)u > +-(p3 - p) 

D 



(3.34) 



For this to hold we need, from ( |3.31| ) 



lim 



L — >oo " 



7 



< u, [e-,-, [cj+i, ejj^2]]u > sin 



37rp 
~2L 



+ < u, (^^/Q{ejejj^i + Cj+iCj ) - 2ej j , u > sin 
= 2p < u, (Lo + Lo)u > +^ (p^ - p) 

D 



Tip 

2L 



(3.35) 



By matching the p and p^ coefficients on the left and right hand side we obtain two 
conditions 



Lim^^oo {3 < u, [cj, [cj+i, ej+2]]u > +a/Q < u, (ejCj+i + ej+iej)u > -2 < u, e^u >^ 
= 

(3.36) 



and 

Linii^^oo (^27 < u, [e^, [ej+i, ej+2]]u > +^/Q < u, {ejej+i + ej+iej)u > -2 < u, eju >j 



-2c7r 



2 / sin_7 ■ ^ 

7 



(3.37) 

Other conditions would be obtained by considering the L dependent corrections to < 
u,u >, which we do not discuss for the moment. From (|3.36|) and (|3.37|) one finds 

TT^c f sino/^"^ 

Li,„,^„ < u, [e,, [e,,„ e,,.]]u [e„ [e,,„ e,,,]U = -- ^— j , (3.38) 

a condition that is in fact independent of u provided u has eigenenergy at distance 0(1/L) 
from the ground state. On the other hand we can have access to some of these averages 
by using the known expression of the free energy of the vertex model (or the Q continuous 
Potts model) in the thermodynamic limit together with the expansion ( p.54|) and the 
relation between the left sides of (|3.37| ) and ( |3.38| ) and conserved quantities. Explicitely 
one has M 



Lim 

where 



^ CX3 /„ \2n-\-l 



CxD 



n=0 

,2n sinh(7r-7)t 



2(2n + 1)! 



In= / r''— H T^dt. (3.40) 

smh Tvt cosh ■yt 



oo 



(3.41) 



So one has in particular 

Coo = sin 7/0 

[cj, [ej+i,ej+2]]oo + cos7(ejej+i +6^+16^)00 = 2sin^7/i 
One finds that QOTI) and ( [CTD imply 

6 _ 247r Jo 48n 

^ " x{x + l) ~ "(x + l)3sin2[7r/(x + l)] {x + 1)^ ^' 
The integrals can easily be evaluated for x integer, and are expressed as rational functions 
of trigonometric functions. One finds that (|3.42[) exactly holds true for x = 1,2,3. For 
X = 3 for instance 

/o = l + -, /i = l + |-. (3.43) 

For X = 4 the righthand side is .701184 instead of .70 while for x = 6 it is .861348 instead 
of .85714. Clearly for x small enough, the result is almost true. 

In the general case, we conclude that the scaling limit of the commutator cannot be 
the commutator of the scaling limits. We need to be more careful to obtain a representation 
of the Virasoro algebra from the lattice models . 



4. Scaling limit and conjecture for the vertex representation 

4.1. Generalities 

The simplest case to study now is the 6 vertex model [Q. In the transfer matrix 
formalism we can think of the various operators so far introduced as acting on the space 
^2L = (C^)^^. The elementary Temperley Lieb matrices Cj, j = 1, . . . , 2L — 1 act as 



e 



+ - - (4.1) 



in the and (j + l)^'^ copies of and as identity otherwise. In the previous equation, 
the basis of is denoted by +, — and E are unit matrices. As explained in [|13| and 
the last Temperley Lieb generator e2L reads generally 

e2L = q-'E+_,+_ + qE_+,_+ - e'^E+_,_+ - e-^^i?_+,+_, (4.2) 

where 99 is an arbitrary complex number. Mathematically the number (p selects different 
representations of the periodic Temperley Lieb algebra (|2.40| ) . Physically it corresponds 
to twisted boundary conditions for the 6 vertex model or the XXZ quantum spin chain in 
the hamiltonian limit. In the following we focus on this limit. Introduce (the normalization 
is such that we now consider the chain as being of length 2L while it was L in the above 
Potts model point of view) the hamiltonian 

2L 

^$^(e,-eoc), (4.3) 

7r sm7 

i=i 

where Coo is ^ independent and ensures that the ground state energy has no extensive 
term. Eq. ( [4.3| ) reads as well 



O /2L-1 _i 



Tc sm7 



J=i . (4.4) 



+ ^^1 ^2L + ^2L + 2 "^^^"^^L - L 2Lec 

Notice that [H^, S^] = where 5"^ = ^cr|. We can therefore decompose 

L 

n2L= nil, (4.5) 



where the sum runs over integers. 
Introduce also P;^ by 

1 / \2 2L 

P^ = -^^— V[e„e,+i]. (4.6) 

The foUowing result has been established by several non rigorous approaches, and carefully 
checked numerically [^[^. Consider the limit L oo. Introduce two other positive 



numbers Tji^Tj such that 

L,Tr,Ti oo, Tr/L tR, Ti/L ti with tR,ti finite. (4.7) 

Then 

lim tr^s^ exp(-TjH^ - iTRp^) = , ] , \^ 

(4.8) 

where the limit is defined in ( ^77|) , e,^ = p = exp(— 7rt/ — iTvtR) and r]{p) = 
pi/24 ]^^^q(1 — p"). This coincides with the trace of 

^Lo-c/24-Lo-c/24^ 

over the space 

0i^a,ao ®^a,ao, (4-9) 
a 

where -Fq,qo is the Gaussian Fock space built of fields of the form 

P{d(j),d^<p,...)expia(l), (4.10) 



where (j) is a chiral boson ||2^ . The sum is taken over e G Z with 

e — e 

= — ^Oi+ + ao + {-)—a-. (4.11) 
The twisted stress-energy tensor reads 

T{z) = ~{d(/}f + iaod^(j), (4.12) 



with 



1 



X + 1 / X 



and we have parametrized q = exp(z7r/(x + 1)). This system has central charge 

c=l-2Aal = l-- -, (4.14) 

and an "effective central charge" equal to one due to the negative dimension operator 
e = 5"^ = 0, = 0, a = ao- 

The same result would hold if instead of (|4.3[) and ([4 .61) we considered any of the 



other (properly scaled) hamiltonians obtained from ( |2.58| ) and (|2.59| ), call them ^ ^ 
and respectively. 

In the following it will be convenient to use the Heisenberg algebra (see eg ) that 
follows from a mode expansion of the chiral boson 

[an, am] = nSn+m, n, me Z. 
Corresponding to the vertex operators expzc/) we have then vectors Va^ao with 



anVoi,oio = for 71 > 0, aot'a,ao = V2(a - 



a,ao 5 



and the Fock space -Fq,qo is built as 



-^a,ao — CcL — m ■ ■ ■ tt—nf^Va^ao- (4-15) 

fc=0 l<Tii <...<nj. 



Then -Fq,q(, is a Virasoro module with generators 

^ oo 

k= — oo 

oo ^ ■ 

Lo = ^ : a-kttk : +-al - 



(4.16) 



k=i 



where normal order is defined by putting the largest index on the right. The space -Fa,ao 
is graded by eigenspaces of Lq 

oo 

^a,ao =0(i^a,ao)„- (4-17) 

n=0 

The dimension of the eigenspace (Fot,oto)^ with Lq eigenvalue — 2aao + n is p{n) the 
number of partitions of n. 



We define the positive definite hermitean form <,> on the Fock space -Fa.ao con- 
travariant with respect to 

(4.18) 



'^n — O'—n-i 



where by convention the vector Va^ao has norm one and Va.,a.o,Va\oLQ are orthogonal for 
a ^ a' . Notice that with this definition the representation of the Virasoro algebra in -Fa,ao 
is not unitary since 



^ oo 



(4.19) 



k= — oo 



Observe that with the natural scalar product in T^fj^, the Temperley-Lieb matrices (^ 
obey also 7^ ej. 



4-2. Scaling limit and Bethe ansatz 

As stressed in section 3 one must be careful with the analysis of 7i|^. First there 
does not seem to be a simple way to imbed Tif^ • Hence it is difficult to follow 

states when L increases and define a suitable limiting process. Also Tif^ is obviously "too 
big" [31]. To proceed further we use the fact that the hamiltonian is diagonalizable 



by Bethe ansatz ||32|| ||33|| . The eigenstates (we restrict to 5"^ > due to the symmetry 
(5"^, 99) {—S^, —'^)) are sums of plane waves with momenta {kj} {j = 1, . . . , L — 5"^) 
solutions of the set of coupled Bethe equations 



2Lkj = 2nlj +ip-^Q{kj,ki), 



(4.20) 



where is the usual kernel of the XXZ chain and {/j} G Z{Z + |) for L — odd 
(even) and are all distinct. The corresponding eigenenergy reads 



E 



7 
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q + q- 



+ cos k. 



2Ler 



(4.21) 



and the momentum obtained by considering properties of the wave function under trans- 
lations reads 



i=i 



(4.22) 



We shall characterize the eigenstates and their corresponding eigenenergies by the asso- 
ciated set of numbers {Ij}- The ground state in every 7i|^ is obtained by choosing the 
maximally packed set of Ij, symmetrically distributed around the origin : the associated 
eigenvalues of H<^ and P,^ reproduce in the limit ( |4.7D the conformal weights of f a,ao with 
a — —^a+ + ctQ + ^q;_. More precisely one has 

L c 
UmL^oo—E{{Ij}) = (a^ - 2aaQ) + (a^ - 2aaQ) - — . 

TT iz 

By shifting as a whole the set {Ij} by e units (to the left or to the right depending on 
the sign of e) one obtains an excited state whose associated eigenvalues reproduce in the 
limit p.7| ) the conformal weights of fa^ao ® Va,ao with a{a) = ^—^aj^ + cto + (— )^a_. 
Creating "holes" in the (shifted) set of integers can lead to further excited states whose 



eigenvalues reproduce in the limit ( |4.7| ) the conformal weights of descendents in -Fa,ao ® 
Fa,aQ- Call {-^j*} the particular distribution that leads to Vct^ao ® Vct^cto- 

CaU {Ij} the 

distribution obtained by making some holes; then the eigenvalues of the associated state 
approaches the ground state as 1/L provided the set {Ij — Ij} has a finite, L independent, 
subset of non vanishing integers . In the simple case when (f = and all the momenta kj 
in ( [4 .201) are real, the corresponding conformal weights read 

2 , (4.23) 

SO the left and right excitations occur for jumps of positive and negative Ij. 

If we characterize the states by the integers {Ij} we can put eigenstates for various L 
in correspondence and define in that fashion a limit process, which we also denote by the 
symbol We can also simply isolate the eigenenergies that will contribute to ( [4. 81 ) . The 
associated set {Ij} has to differ from a symmetric, maximally packed distribution by finite 
global shifts and finite numbers of holes, these deviations being moreover L independent 
35[] . This picture becomes transparent in the case q = i where the O term in ( [4.20| ) 



vanishes and the spectrum is the one of a free theory. See also the Ising model of the 
previous section. We sometimes call states characterized by such distributions of integers: 
scaling states. We shall talk also about the limit of these states, meaning that we consider 
the family of states characterized by similar patterns of {Ij} as L increases. Usually the L 
dependence will not be written explicitely, in particular in the numerical results of section 
5. 



Note that we implicitely assumed that aU eigenenergies could be found from the Bethe 
ansatz. This is generaUy beheved to be true. For special values of q or (f> this may be true 
up to degeneracies. But in such cases, using the known additional symmetries of the 



hamiltonian , the correct multiplicities can be recovered. 

Note that in practice, when the hamiltonian is diagonalized numerically, it is easy to 
"follow" the eigenenergies and their associated eigenstates by simple order. Comparing 
the k^^ {k fixed) eigenenergies for increasing values of L, is a reliable limit process. Notice 
also that the momentum can be measured without finite size corrections from V, which 
makes the identification of terms in (|4.8|) easier. 



4.3. Lattice Virasoro generators 

We now state the main conjecture of this section. Introduce the lattice quantities 

= 7^ I ^— E ""'^'^"^ f^^- - + ^^[e.' e.+i]) i + 7^<Jn,o, (4.24) 

2ti 7rsm7^ \ 7rsm7 y 24 



and 



27r 1 7rsm7^ \ 7rsm7 / J 24 

Consider the following double limit process. Choose a value of a and a and a pair of 
integers A^, N. Using the characterization of Bethe eigenstates by integers, or simply by 
ordering the eigenenergies, select for every L the set of 

N N 
n=0 n=0 

eigenstates associated with gaps gaps which, after multiplication by Tr reproduce in the 
limit L — » 00 the conformal weights /ia.ao + ^ ^ind /iq^qo + n with n < N,n < N (for L 
large enough but finite all these states will be present: we implicitely suppose we always 
are in such a situation). Consider the action of the generators /„ (resp. In) restricted to 
these states. We conjecture that it furnishes in the limit L — 00 a representation of Vir 
(resp. Vir) restricted to 

N N 

^ ® (-^a,ao)n ® {F(i,ao)fi ' 
n=0 n=0 

that is more precisely a representation of PVirP (resp. PVirP where P is the projector 
on the above subspace). Moreover we conjecture that the natural scalar product in TC2L 
(with < +, + >=< — , — >= 1 and < +, — >=< — , + >= 0) coincides in the L ^ 00 limit 
with the hermitean form defined previously in terms of Heisenberg algebra. 



5. Numerical and analytic checks 

A first remark is that 

[ln,V] =nln, ^n,V] 



-nL 



(5.1) 



so the lattice InJn have non vanishing matrix elements only between eigenspaces of the 
lattice momentum V whose eigenvalues differ by in. This is of course the same result as 
what is expected in the L ^ oo limit. 

5.1. Level one 

To start consider the ground state in each sector H^l; call it u. It should correspond to 
Va,ao with a = —^a+ + ao + ^a-. In each sector Tif^ consider also the first (normalized) 
excited state with momentum +1, which should correspond, up to a proportionnality 
factor, to a_ifa,ao; call it v. In the continuum theory one has 



and thus 



e 



We thus checked numerically in a variety of cases that 



"0 - ^ ~2^~ 



and similarly 



We also checked that 



(v,/_iu) ^ V2 
(u, hv) -V2 



an a+ H a_ 

"2+2 



(v,riu) ^0, (u,Liv) ^0. 



A special value is = —2aoa- and = for which 

(v,Z_iu) -^0, 
and = 2aoCK- and = for which ( |5.5D becomes 

(u,/iv) ^0. 



becomes 



(5.2) 
(5.3) 

(5.4) 

(5.5) 
(5.6) 
(5.7) 
(5.8) 



In fact the first result holds even for finite L. It is a consequence of the representation 



theory of the periodic Temperley Lieb algebra this is explained in subsection 5.7. 
Some numerical results for the foregoing quantities are given in tables 1,2,3,4. 
Recall finally that from (157^) we have exactly 

(v,r±,u) = 0, (u,Z±„v) = 0, (5.9) 

for n > 1. 



5.2. Level two 

Consider for simplicity the case Tifi^*^ only, and identify the two orthonormal excited 
states above the ground state u, with momentum two. For any x it is easy to find one 
of them by numerical solution of the Bethe ansatz equations. It corresponds to integers 
{Ij} obtained from {Ij} by shifting the two rightmost integers by one unit. We call it wi. 
We had difficulties finding the other one by numerical solution of ( |4.2U| ) (which should be 
obtained by shifting the the rightmost integer by two units to the right) so we used direct 
diagonalization (Lanczos algorithm |]H^). Call this other state W2. In the large L limit we 
expect 

|wi , „ = (xa^i +ya-2) ka,ao >, (5-10) 

a/2x^ + 2y^ 

where x, y are coefficients to be determined. To help us determine the coefficients y we 
made the following observations. Consider first the the case e^p = 0. We found then 

< U, /2W1 >= - < Wi, /_2U >, 

this result being true for any finite L, and a fortiori for L 00. From this we deduce 

a: = Oife<^ = 0. (5.11) 



Hence 



and 



|wi >t^ —a-2\Va,ao K2 >^ --j=a_i\Va,o,o (5-12) 

< wi, /_2U >^ 4ao 

1 . (5.13) 



We also observed that and 



< W2, /-2U >- 

v2 



< U, /2W2 > = < W2, I-2U >, 



holds exactly for finite L too. Corresponding numerical results are shown in table 5. 
Similarly in the case = —2aoa- we found 

< Wi, /_2U >= 0, 

this result being true for finite L. From this we deduce 



X = —2V2aoy if e^p = —2aoa-, 



(5.14) 



and then 

|wi >t-^ , = (-2A/2aoa-i + a-2) \va,ao >, (5.15) 

A /2 + 16aQ ^ ^ 



and 



|W2 >^ , = fa^l + 2^2000-2) \Va,ao >= ^^===^-^-2^0,00 > • (5.16) 



So we get the results 

, 8a, 
< u, /2W1 >^ - 



< u, /2W2 ^===% . (5.17) 

a/2 + 16a5 

< W2, /-2U >^ y — ^ — 

Corresponding numerical results are shown in table 6. 

We can also consider action of l2- An especially interesting quantity is then (still for 
— — 2aoa-) 

< U, ^2^1 >< Wi, /_2U > + < U, /2W2 >< W2, l-2n >^ -. (5.18) 

Numerical results are given in table 7,8. They actually provided very good estimates of 
the central charge. 

We can also consider action of l'^^. For instance, for = 0, 

< wi, /_iv >< V, /_iu 2ao 

< W2, f-iv >< V, l-iw >^ 2v2aQ 

while for = —2aoa-, 



< u, /iv >< V, /iw >^ 4q;oY 1 + Sag- (5.20) 

Corresponding numerical results are in table 9. 

Also, at level two we can also consider action of l-i and l-i. For instance consider 
the eigenstate |w > with vanishing momentum, associated with a gap that reproduces the 
conformal weights h = h = 1 in the limit ( [1.8|) . One expects then 

< w, Liv >< V, /_iu >— > 2aQ. (5.21) 



Numerical results for this quantity are given in table 10. 

At level two we can finally test numerically a non-trivial case of degenerescence. Con- 
sider, with e^p = —2aoa- and = the state |u' > of conformal weight h2i. This state 
is degenerate at level two. From the Coulomb gas mapping we expect that 

< w', /_2u' > -—-A—- < w', Uy' >< v', /_iu' >, (5.22) 
Z[zn2i + J- j 

will vanish for large systems, where |v' > and |w' > are the appropriate excited states. 
Numerical results are shown in table 11 and converge very well to zero. 

5.3. Level three 

Consider now Tif^^*^ and identify the three orthonormal excited eigenstates states 
above the ground state, with momentum three. Call them yi,y25y3- In the large L limit 
we expect 

|yi > ^ . ^ ,^ =^(C1Q-1 + c[a-ia-2 + a-3)\Va,ao > 

V6cf + 2c^^ + 3 

1^2 > ^ 2 \ r ^('^2a-i + a_ia_2 + C2a_3)|wa,ao > (5.23) 

Y DC2 + 2 + 0C2 

1^3 > ^ /a 2 \ ^2 ^*^^30-1 + C3a_ia_2 + a-3)\Va,ao > ■ 

Y 0C3 + 2C3 -|- o 

To determine the unknown coefficients, we use orthogonality of the states, which implies 

6ciC2 + 2ci + 3c2 = 

6C3C2 + 2c3 + 3c2 = (5.24) 

6C1C3 + 2ciC^ + 3 = 0, 

and the numerical obervation that for = 0, 

< yi, /-3W > = - < u, hy-^ > 

< 72,1-3^ > =< u, ^372 > (5.25) 

< 73, 1-3^ > = - < u, Z3y3 > . 

From the above, we deduce that all the coefficients vanish except for ci , C3 which are related 

by 

C1C3 = -1/2. (5.26) 



In addition, the above matrix elements of /_3 are now given by 

<y2,/_3U> = (5.27) 
6-\/3ciao 

There does not seem to be any way to determine the value of the unknown c\. That 
C1C3 7^ shows that the eigenstates of finite chains do not converge in general to pure 
monomials in the Heisenberg algebra as was the case at level two. Even if c\ is not known 
observe from ( ^.27|) that 

(< yi, /-3U >2 + < y3, Z_3U >2) = 3v^ao, (5.28) 

a result in dependent of c\ that can therefore lead to numerical study. Numerical results 
are shown in table 12. 

Similarly for e^p = —2aoa-, one of the three excited states y2 corresponds, in the 
continuum limit, to the level 3 descendant of the vacuum |t'a,ao > and is given by 

^ ^ ° (5.29) 

= , =(a-ia-2 + 2A/2aoa-3)ka,ao > • 

V 2 + 24^2 

While the other two excited states yi, 1/3 correspond, in the continuum limit, to two level 
two descendants of the null states and are therefore given by 

\jl >^ I (^2 + -^1^1 )«-l ka.ao > 

^(3 + 54a2)(i + 2Ai)2 + 6(i + 8a2)2 
1 



(ciai?_i - 3\/2aoa-ia_2 +a_3)|'i;a,ao > 

^3 + 36aQ + 6cf 

|y3 >^ I (4 + ^3l?)a-i\va,c.o > 

^(3 + 54a2)(i + 2A3)2 + 6(i + 8a2)2 

> 



(5.30) 



v/3 + 36al + 6cl 

where we have traded the unknown Ai(3) by another ci(3) for convenience. Orthogonality 
then implies that 

C1C3 = -(l/2 + 6a2). (5.31) 



We therefore expect 

< yiJ-3U > — > 
< kVi > — ^ - 



3 + 36a^ + 6cf 



< y2j-3U > — > \/2 + 24a^ 
< u, hy2 > — 



2-48a2 , (5.32) 



2 + 24^2 
< ysJ-su > — >0 

ISv^ao 

< u, hy3 > — ^ 



v/3 + 36^2 + Qcj 

Numerical calculation shows that < yi(3),/_3W > vanish exactly for finite L and as in the 
case of e^p = 0, the coefficient ci remains unknown. However 

(< y„ L3U >^ + < y3. i-3U = -5%,, (5.33) 

a combination that is independent of ci. Some numerical values are given in table 12 and 
13. 

5.4- Higher order approximations 

According to eqs. ( ^.581 ) and ( |2.59[ ) we can form lattice approximations to the Virasoro 
generators by choosing any generic term of higher order conserved quantities. As an 
example we give in tables 14 and 15 measures of the central charge from ( ^.18|) . The 
results are comparable to the ones obtained with the lowest hamiltonians. 

5. 5. The double limit process 

The double limit process is very important in the conjecture. However it plays only a 
small role numerically, as was already observed in the discussion after eq. ( p.42|) . Consider 
for instance the lattice quantity 

< /_iu, /_iu > . (5.34) 



Compare this expression to ( |5T^ ) . For instance for e<^ = and = eq. ( |5T^ ) gives 
2aQ = 2x{x+i) ■ ■^'-'^ ^'^ ~ ^2aoCi- eq. (|5.2| ) gives zero. Measures of (|5.34| ) for x = 3, 7 are 



given in tables 16 and 17. The data do not behave as nicely as the ones for 



< /_iU,V >< V,Z_iU >, 



(5.35) 



which can even give exact results (as a consequence of Temperley-Lieb representation 
theory) for finite L. It is very hkely that they are not converging to ( p.2| ) since they are 
monotonous and above ( |5.2D for L > 7. This is expected for the reason aheady discussed 
in section 3.2: in ( |5.34| ) aU intermediate states contribute, including those at large distance 
from the ground state. Although matrix elements of /_i between u and states different 
from V are expected to vanish in the limit L oo, the sum of all the corresponding 
contributions may well remain finite. 

5.6. Numerical Comments 

In general, we think that the numerical agreement with the conjecture is very good. 
The measure of c from ( |5.18D in particular provides results of the same precision as direct 
study of finite size effects for the ground state. It so happens that for a particular quantity 
and a particular value of x the agreement is not perfect (see eg table 1). This however can 
usually be attributed to data being non monotonic, hence difficult to extrapolate. But for 
a given value of x, most tests give good results. 

The difference between limit of commutator and commutator of limits seems obeserv- 
able numerically, but is rather weak. 

5. 7. Some analytic checks from Temperley-Lieb representation theory 

Some of the conjectured results actually hold exaclty for finite L and can be established 
by using the representation theory of the periodic Temperley-Lieb algebra, in particular 
the analysis of the vertex model representation in |2^. Consider first the case x irrational. 



Then for generic the space provides an irreducible representation of the periodic 
Temperley-Lieb algebra. This representation however breaks when is a multiple of 27 
(recall 7 = or equivalently e^^ is a multiple of —2aQa-. Setting 

= 2n7, n > ^^ (5.36) 

Tif^ has an irreducible component -Rf^'" of dimension {^2l ) ~ (^21") • "^^^ generating 
function (|4.7|) reduced to this representation reads 



Lim exp(-T7H^ - iTRp^) = T^' - Tq. = (5.37) 



where recall that JF^ refers to the right hand side of (|4.8|) . Eq. ( ^.37|) corresponds to the 
trace of 

^Lo-c/2^-Lo-c/2A^ 



over the Vir ® Vir irreducible component of (O). The matrix elements of the lattice 



_______ _______ QZ ^ 

generators ( |4.24D and ( |4.25D between states in the irreducible component -R2L' states 
out of it vanish. As an example, the matrix elements of the lattice Virasoro generators 
vanish exactly between the ground state u (that belongs to the irreducible component) 
and the excited states v and wi when 5"^ = and e^p = — 2cKoct- (that do not belong to 
it) . This is an exact lattice analog of the cancellations in the action of Virasoro generators 
deduced from ( |4.12| ) in Fq, q,q. 



In the case when x is rational , the representation theory of the periodic Temper ley- 
Lieb algebra becomes more involved, and the space -R2L' becomes itself reducible. Restrict 
for simplicity to the case x + 1 integer. A new irreducible component P2L^ appears then 
of dimension 

^ /L- + l{x + l)\ /L-n + Z(x + l 
2L )'[ 2L 

where the sum truncates for negative arguments in the binomial coeffcients. The generating 
function p.7| ) reduced to this representation reads 



x-l 

Lim tr^s-,r^ exp(-T7H^ - iTRp^) = ^ Xr,n-s- Xr,n+s- , (5.38) 

r=l 



and corresponds of course to the trace of 



^Lo-c/24-Lo-c/24^ 



over the irreducible component of (^]^) that is 

x-l 

^ Vir^,„_szVir^,n+s-, (5.39) 

where in ( |5.38| ) Xrs denotes the character of the Virasoro algebra in the irreducible repre- 
sentation Vir^-s and labels r, s are Kac labels parametrizing the conformal weights 

[(x + l)r — xs]'^ — 1 



h 



rs 



4:X{X + 1) 

see and references therein. 

Matrix elements of Temperley-Lieb generators, hence in particular of the lattice Vi- 
rasoro generators, vanish then exactly between states in pfj^'" and states out of it, in 
agreement with corresponding results for Virasoro generators in -Fq, q.^. To analyze fur- 
ther the irreducible representations in the rational case it is better to turn to the RSOS 
representation. 



6. Scaling, conjecture and numerical checks for the RSOS representation 
6.1. Generalities 



We restrict for simplicity to x integer. We recall the RSOS representation of the 
Temperley-Lieb algebra. The space H^l^^ is now spanned by vectors . . . , I2L > where 
the li take values = 1, . . . , a; with the constraint — \ = 1. This space has dimension 
Trc^^ where c is the incidence matrix of the An,. Dynkin diagram. One has 

=^ih, I'l) ■ ■ ■ l'j_i)d{lj-i, lj+i)d{lj+i, I'j^^) . . . d{l2L, 1'2l) 

(sin(7r/,/(x + l))sin(7r/^./(x + l))'/^ ' (6-1) 



sin(7r/j_i/(x + 1)) 



with I2L+1 = h- Define 



and 



2L 

H = ^ V(e, -eoo), (6.2) 

TT sm 7 

P = -f^^) i;[e„e,+i]. (6.3) 
The space 7^^^*^^ decomposes as 



s=l 



where the irreducible representations of the periodic Temperley-Lieb algebra. 

One has from subsection 5.7 



x-l 



Lim tr Q,s exp(-T/H - iTrP) = XrsXrs- (6.5) 

r=l 

With the natural scalar product in Ti^f^'^^ one has now = tj. 
6.2. Scaling limit 

It is well known that the eigenstates of (|6.2| ) are a subset of the eigenstates of (^ 



for appropriate choices of the spin 5"^ and the twist angle The same scaling limit 

can therefore be defined as in the vertex case. 



6. 3. Conjecture 

We now state the main conjecture of this section. Introduce the same lattice quantities 



[ej,ej+i] ) > + 7^5n,o, (6.6) 




7rsin7' ^' ' J 24 



. V + -^5„,o. (6.7) 

7rsm7 / 24 

Consider the foUowing double limit process. Choose a value of s and a pair of integers 
iV, N . Using the characterization of Bethe eigenstates by integers, or simply by ordering the 
eigenenergies, select for every L the set of gaps which, after multiplication by reproduce 
in the limit L ^ oo the conformal weights /ia,ao + ^ and /iq^qo + n with n < N,n < N (for 
L large enough but finite all these states will be present: we implicitely suppose we always 
are in such a situation). Consider the action of the generators 1^ (resp. /„) restricted to 
these states. We conjecture that it furnishes in the limit L — oo a representation of Vir 
(resp. Vir) restricted to 



TV N x-1 



n=0 71=0 r=l 

that is more precisely a representation of PVirP (resp. PVirP where P is the projector 
on the above subspace). Moreover we conjecture that the natural scalar product in Ti^f 
coincides in the L —>■ oo limit with the hermitean form for which Z+ = 1-^. 

6.4- Numerical checks 

Numerical checks are much of the same nature as the ones given for the vertex case. 
We shall give just a few examples. 

First consider the state at vanishing momentum whose gap reproduces in the limit 
( |675| ) the weights h = 2,h = 0. Call it |w2 >. One has thus 



|w2>-,/^-^L_2|u>, (6.8) 

where by |u > we denote the ground state of the RSOS model. By Virasoro commutation 
relations one expects 

< W2,/-2U >^ (6.9) 
Results for this quantity are given in table 18. 



7. The case of fixed boundary conditions 

So far we have deah with systems without boundaries. This is the most favorable for 
numerical checks. On the other hand we can only select combinations of left and right 
Virasoro representations of the type 

x-l 



^ Vir^sVir^s/, (7.1) 



r=l 



using the lattice symmetries. The irreducible representations of left or right algebras have 
to be selected by hand after the hamiltonian has been diagonalized, which is not too 
satisfactory. This problem can be solved by turning to fixed boundary conditions. Indeed, 
following consider the RSOS model of section 6 restricted to the space 'H'^l'^ spanned 



as before by vectors . . . , I2L > with 

h = a, I2L-1 = b, I2L = c. (7.2) 

Then one has 

Lim tr a/bc exp(-TH) = Xda, (7.3) 

''■2L 

where 

d = mf{b,c). (7.4) 

To get all the posible parities of character labels one can also consider the RSOS model with 
an odd number of heights where the same formula holds. We can therefore easily select 
in the lattice model a space whose scaling states will correspond to a single representation 
of the Virasoro algebra. The complete Virasoro action can be conjectured based on the 



arguments of section 2 using the same lattice stress energy tensor and the results of ||4C1|| . 
One expects that 

( 2L-1 . , x2 2L-2 . 1 

2L I 7 / N , / 7 \ r 1 . '^JTT I c 

^- = -^^^1^^ - ^ + [V;^ ) 1. ^ j + 24'^-°' 

(7.5) 

tend to the Virasoro generators with the limit process already outlined in sections 4 and 
6. We have preformed detailed numerical checks of this conjecture as before. The conver- 
gence is not as good due to the free boundary conditions but the results are satisfactory. 
Examples are given in table 19 to be compared with table 18. 



8. Conclusion 



The double limit process in the scaling can probably be overcome by considering 
operators at different "times" . Indeed consider for instance 

UT) = exp (^TT^n^ e-^^lne^^, (8.1) 

for say the RSOS model. Consider some scaling state |u > and act on it with 

ln{Ti)l^{T2), (8.2) 

with Ti < T2 and consider the limit Ti,T2, L — > 00 with T^2)/L — > ii(2)- We expect that 
acting on |u > will give the following results. The first term e^^H produces simply a factor 
g7rT2/L(/i„+/i„)|^ ^ before action of 1^ couples |u > to many other states. Except for 
the state that corresponds to L^|u > all the matrix elements vanish for large L, but before 
their sum could give finite contribution after insertion of Ln- However we now have the 
exponential factor e(^i~^2)H -por non-scaling states the value of H is finite and therefore 
in the limit Ti,T2 ^ 00 their contribution is damped out exponentially fast. For scaling 
states the scries with generic term converges and therefore limit can be taken 

term by term. So we expect that all the unwanted terms will disappear to give simply the 
state corresponding to L^L^ |u >. Hence the reformulation 

LimTi,T2,L^oo {ln{Ti)lm{T2) - /^(Ti)/„(T2)) |u > = 

/ c o \ • (8.3) 

Limy^i^oo ({ri - m)ln+m{T) ~ ri)Sn+m] |u > 

The approach we have used is simple minded but we think it gives nice results. We 
have observed with remarkable precision the built-up of the central term , the null vectors 
structure, the metric properties, and we believe have indeed a way of extracting the Vi- 
rasoro algebra from the lattice model. Numerically, the integrable lattice models provide 
probably the best regularization for quantum field theories in 1-1-1 dimensions, and we are 
not aware of works using more direct regularizations (eg. discretization of a free boson) 
that would give comparable results. To make this more interesting analytically one would 
need to reproduce our computations exactly using the Bethe ansatz. This seems difficult 
but maybe not impossible. If the program can be carried out, it will provide a bridge 
between the Virasoro algebra and lattice integrability. In mathematical terms, we have 
a conjectured homomorphism between the universal enveloping algebra of the Virasoro 



algebra and the Temperley-Lieb algebra when the number of generators goes to infinity, a 
very intriguing algebraic result. 
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Figure captions 



Figure 1: Conventions for labelling the sites and edges of the square lattice and its dual 

Figure 2: A Potts model on the square lattice with ^ Ky can be considered as dis- 
cretization of an isotropic continuum medium by rectangles. We characterize the rectangles 
by the anisotropy angle 9. 

Figure 3: The operator td propagates in the direction of the light arrow. The operator 
Vtd propagates in the direction of the thick arrow. 

Figure 4: When ^ Ky, Vtd depends both on Lq + Lq and Lq — Lq. 



Table 1 Numerical values of — < u, hv > with 6,^ = in the = sector." 



2L 


X = 3 


4 


5 


6 


7 


8 


0.20266 


0.15784 


0.12955 


0.11000 


0.094502 


10 


0.20330 


0.15818 


0.12974 


0.11011 


0.095719 




U.zUoDi 


U.ioooz 


u.izy /y 


nil ni o 
U.llUlz 


u.uyo/Ui 


14 


0.20379 


0.15838 


0.12979 


0.11009 


0.095656 


16 


0.20390 


0.15841 


0.12977 


0.11005 


0.095606 


extrapolation 


0.20416(2) 


0.1584(3) 


0.1292(2) 


0.1092(3) 


0.0945(3) 


conjectured 


0.20412 


0.15811 


0.12910 


0.10911 


0.094491 


Numerical values of < v, Z_iu > 


are identical. 


Conjectured numbers 


are V^CKq 




Table 2 Numerical values of 


— < u, liv > with = —2a 


oCK- in the 


= sector.'' 


2L 


X = 3 


4 


5 


6 


7 


8 


0.40289 


0.31422 


0.25817 


0.21938 


0.19088 


10 


0.40475 


0.31521 


0.25874 


0.21973 


0.19109 


12 


0.40577 


0.31570 


0.25897 


U.ziyoz 


U.iyiii 


14 


0.40639 


0.31597 


0.25906 


0.21982 


0.19106 


16 


0.40679 


0.31612 


0.25909 


0.21978 


0.19098 


extrapolation 


0.4082(2) 


0.31638(5) 


0.2591(2) 


0.2183(3) 


0.1899(3) 


conjectured 


0.40825 


0.31623 


0.25820 


0.21822 


0.18898 



< V, Z_iu > is exactly zero by TL representation theory. 



Table 3 Numerical values of < v, Z_iu > in the = 1 sector.*^ 



2L 




= 




2q!oQ;_ 


X = 3 


7 


X = 3 


7 


8 


0.37781 


0.51645 


0.55468 


0.60132 


10 


0.38728 


0.52838 


0.57403 


0.61729 


12 


39282 


53742 


58494 


62681 


14 


0.39634 


0.54066 


0.59170 


0.63299 


16 


0.39874 


0.54416 


0.59619 


0.63729 


extrapolation 


0.4082(2) 


0.547(2) 


0.6123(4) 


0.660(4) 


conjectured 


0.40825 


0.56695 


0.61237 


0.66144 



The conjectured value is V 2|q!o — -^OL^ + 



Table 4 Numerical values of — < u, l\\ > in the = 1 sector. 



2L 




= 




—2aoa- 


X = 3 


7 


X = 3 


7 


8 


0.79920 


0.72006 


0.97543 


0.80498 


10 


0.80561 


0.72945 


0.99188 


0.81843 


12 


0.80922 


0.73516 


1.00097 


0.82630 


14 


0.81140 


0.73901 


1.00646 


0.83130 


16 


0.81281 


0.74150 


1.01001 


0.83469 


extrapolation 


0.8169(3) 


0.74632(1) 


1.021(1) 


0.849(2) 


conjectured 


0.81650 


0.75593 


1.02062 


0.85042 



The conjectured value is V 2|ci;o ~l~ ^CK-)- — 



Table 5 Numerical values of < tui, I-2U > and < W2, 1-2U > with = 0.® 



2L 


< W2i I 


-2W > 




-2U > 


x = 3 


7 


a; = 3 


7 


8 


0.70117 


0.70264 


0.55799 


0.26624 


10 


0.70357 


0.70491 


0.56537 


0.26839 


1 9 

±Zr 


u. 1 u^ou 


U. 1 UUIZ 




n 9fiQ49 


14 


705fi2 


70680 


571 68 


26QQ2 


16 


0.70608 


0.70721 


0.57318 


0.27017 


extrapolation 


0.70705(2) 


0.7073(3) 


0.57753(3) 


0.2704(2) 


conjectured 


0.70711 


0.70711 


0.57735 


0.26726 


" In the =0 


sector, <^ u, ^2*2 


> = < W2, /-2U 


>, and < u, /2WI 


>= — < wi, I 



Cojnectured values are respectively —7^ and 4q;o. 



Table 6 Numerical values of < W2, i-2U >, < u, Z2W2 > and — < u, Z2W1 > 
with e<p = —Iolqol-J 



2L 


< W2, 1-2^ > 


< U, Z2W2 > 


- < U, Z2W1 > 


a; = 3 7 


a; = 3 7 


a; = 3 7 


8 

10 
12 
14 
16 


0.76007 0.71725 
0.76243 0.71922 
0.76357 0.72022 
0.76415 0.72075 
0.76445 0.72104 


0.31235 0.60688 
0.31725 0.61099 
0.32000 0.61331 
0.32171 0.61475 
0.32285 0.61571 


1.01588 0.51990 
1.03427 0.52490 
1.04453 0.52738 
1.05082 0.52869 
1.05495 0.52940 


extrapolation 


0.76479(1) 0.72140(1) 


0.3274(6) 0.6200(7) 


1.0692(3) 0.53029(2) 


conjectured 


0.76376 0.71962 


0.32733 0.62037 


1.06905 0.52523 



In the — sector, < wi, i_2u > — exactly. The conjectured results are respectively 



l+8ag l-24ag Sap 
2 ' ^2+16a2' ^l+8a2' 



Table 7 Numerical values of < u, /2W1 >< wi, Z_2U > + < u, Z2W2 >< W2, i-2U > 
with = and = O.^ 



2L 


a; = 3 


7 


8 


0.18029 


0.42282 


10 


0.17537 


0.42487 


12 


0.17269 


0.42601 


14 


0.17107 


0.42670 


16 


0.17002 


0.12715 


extrapolation 


0.1667(3) 


0.4288(2) 


conjectured 


0.16667 


0.42857 



The conjectured value is 1/2 — 16ci!q. 



Table 8 Numerical values of < u, l2^i >< wi, Z_2U > + < u, Z2W2 >< W2, 1-2^ > 
in the = {) sector and e^p = —2aQa-.^ 



2L 


a; = 3 


4 


5 


6 


7 


8 


0.23741 


0.33642 


0.38711 


0.41659 


0.43529 


10 


0.24188 


0.34123 


0.39172 


0.42096 


0.43944 


12 


0.24434 


0.34388 


0.39426 


0.42336 


0.44172 


14 


0.24583 


0.34550 


0.39580 


0.42480 


0.44308 


16 


0.24680 


0.34655 


0.39680 


0.42573 


0.44395 


extrapolation 


0.2501(3) 


0.3501(3) 


0.4002(3) 


0.4287(3) 


0.4467(2) 


conjectured 


0.25000 


0.35000 


0.40000 


0.42857 


0.44643 



The conjectured value of < u, Z2WI X wi, i_2u > + < u, Z2W2 X W25 1—2^ ^ is equal to C 



Table 9 Numerical values of < u, liv >< v, ZiWi > and < u, liv >< v, I1W2 > 



21 


< u, Ziv X V, Ziwi > 


< u, Ziv >< V, Z1W2 > 


e<^ = 




e<p = 


X = 3 7 


X = 3 7 


X = 3 7 


8 

10 
1 9 

±Zr 

14 
16 


-0.24115 -0.11157 
-0 25788 -0 120Q5 
-0.26721 -0.12535 
-0.27291 -0.12798 
-0.27663 -0.13188 


0.51308 0.22911 
55112 24558 
0.57256 0.25464 
0.58576 0.26009 
0.59444 0.26359 


0.059246 0.013281 
0.059256 0.013351 

0.059229 0.013231 
0.059208 0.013188 


extrapolation 


-0.289(1) -0.1303(1) 


0.625(3) 0.2714(2) 


0.05895(3) 0.01282(3) 


conjectured 


-0.28868 -0.13363 


0.62361 0.27199 


0.058926 0.012627 



* For e^p = 0, numerical calculation shows < wi, l — iv X v, Z_iu — ^ ^1'^ ^iwi > 

and < W2, l—iv X v, /_iu >=< u, Ziv X v, ZiW2 >• While for = — IcHqOI—, 
< wi(2)) l—l'v X V, Z_iu > and < u, Ziv X v, liw2 > vanish exactly. This is consistent 

with |v > > Z/_i IVq^qq > and |w2 > ^ 2 1 '^01,0:0 ^ respectively. The conjectured values 

are respectively —2ao, Aaoy^l + SO-q, 2-\/2Q!q. 



Table 10 Numerical values of < w, /_iv X v, /_iu > with Cip — and — 0.^ 



2L 


X = 3 


7 


8 


0.040467 


8.9975E-3 


10 


0.040377 


8.9183E-3 


12 


0.040442 


8.8122E-3 


14 


0.040540 


8.7425E-3 


16 


0.040639 


8.6926E-3 


extrapolation 


0.04157(2) 


8.713(7)E-3 


conjectured 


0.041667 


8.9286E-3 



Conjectured value is 2cI!q. 



Table 11 Numerical values of < w',Z_2u' > ~ 2(2/^21+1) ^ ^-1'^' X ^-i^' > 
with = — 2q;oq;-. 



2L 




a; = 3 




7 


8 




0.25098 




0.21854 


10 




0.17640 




0.14104 


12 




0.12864 




0.09025 


14 




0.09722 




0.05808 


16 




0.07576 




0.03710 


extrapolation 




-0.00031(2) 




-0.00047(3) 


conjectured 












1/2 

Table 12 Numerical values of (< u, Isy^ >2 + < u, /ayg >2) ^ . 






= 




= — 2Q!oQ!_ 




2L 


X = 3 


7 


a; = 3 


7 




8 


0.98806 


0.48071 


1.71282 


0.92484 




10 


1.01080 


0.48505 


1.76984 


0.93693 




12 


1.02535 


0.48884 


1.80542 


0.94619 




14 


1.03469 


0.49135 


1.82827 


0.95221 




16 


1.04092 


0.49297 


1.84364 


0.95614 




extrapolation 


1.0512(5) 


0.4958(1) 


1.852(8) 


0.9634(1) 




conjectured 


1.06066 


0.49099 


1.89737 


0.95669 





k Conjectured values are respectively 3\/Qoiq and — y===^=^ for e<n — and — — 2(Xq(X—. 



Table 13 Nmnerical \-alucs of < u, /:jy2 > and < /-.-jU > . 



2L 


< u, l3y2 > 


< y2>^-3U > 


= 




= — 2q!oQ;_ 


x = 3 7 


a; = 3 7 


x = S 7 


8 

10 

1 

14 

16 


1.32379 1.32759 

1.00o4/ l.oODOD 

1.37152 1.37431 
1.38284 1.38543 

1.39030 1.39275 


0.53730 1.11815 

0.58933 1.17635 
0.60025 1.18940 

0.60710 1.19803 


1.50860 1.37583 
1.53200 1.41013 

1 c:/l7'?'? 1 /1 17/1 /I 

1.55707 1.42766 

1.56316 1.13139 


extrapolation 


1.3986(7) 1.4016(4) 


0.6314(2) 1.204(2) 


1.57505(1) 1.457(1) 


conjectured 


1.41421 1.41421 


0.63246 1.23017 


1.58114 1.45160 



Numerical calculation shows for Co:, — < y2, / — 3u > — < u, lsy2 Conjectured values for the 



above are respectively V2, 7-^°"° and V2 + 24a^. 



Table 14 Numerical values of < u, /2W1 >< wi, I-2VL > + < u, /2W2 >< W2, /-2U > 
computed using h^'^\ in the — sector and e^p = — 2q;oQ;_. 



2L 


X = 3 


4 


5 


6 


7 


8 


0.20264 


0.28469 


0.32633 


0.35048 


0.36577 


10 


0.21879 


0.30716 


0.35190 


0.37780 


0.39418 


12 


0.22797 


0.31987 


0.36630 


0.39313 


0.41008 


14 


0.23366 


0.32770 


0.37514 


0.40252 


0.41980 


16 


0.23741 


0.33286 


0.38094 


0.40867 


0.42511 


extrapolation 


0.2.50(1) 


0.351(3) 


0.101(2) 


0.129(2) 


0.1332(1) 


conjectured 


0.25000 


0.35000 


0.40000 


0.42857 


0.44643 



Table 15 Numerical values of < u, /2W1 >< wi, I-2VL > + < u, I2W2 >< W2, i-2U > 
computed using h^^^ in the = sector and e^p = —2aoa-. 



2L 


X = 3 


4 


5 


6 


7 


8 


0.17297 


0.24065 


0.27457 


0.29412 


0.30645 


1 n 


u.iy / oy 


U.Z / D14 


U.01041Z 


U.oooUO 


U.OOZoZ 


12 


0.21270 


0.29718 


0.33959 


0.36400 


0.37938 


14 


0.22209 


0.31050 


0.35486 


0.38038 


0.39645 


16 


0.22837 


0.31940 


0.36507 


0.39131 


0.40783 


extrapolation 


0.252(2) 


0.353(3) 


0.403(4) 


0.432(4) 


0.449(4) 


conjectured 


0.25000 


0.35000 


0.40000 


0.42857 


0.44643 



Table 16 Comparision of < /_iu, Z_iu > and < /_iu, v >< v, Z_iu > with a; = 3 
and = 0."^ 



2L 


e<p = 






< /_iU, V >< V, /_iU > 


< /_lU, /_iU > 


8 


0.041248 


0.041069 


6.0139E-5 


10 


0.041579 


0.041330 


1.0574E-4 


12 


0.041744 


0.041459 


1.3536E-4 


14 


0.041833 


0.041531 


1.5382E-4 


16 


0.041882 


0.041575 


1.6498E-4 


extrapolation 


0.041946(2) 


0.041656(2) 


1.800(l)E-4 


conjectured 


*** 


0.041667 


*** 



The numerical value of < lu, v X v, Z_iu > is exactly zero for 6^, — — loioOi—. For e^, — 0, 
the conjectured value is 2q!q. 



Table 17 Comparision of < I 


_iu, Z_iu 


> and < Z_iu, 


V >< V, Z_iu 


> with X = 7 


and = 0. 























— Soi^Qi 




2L 


< l-lU, l-iU > < L 


-lU, V >< V, l- 


lU > 


< l-lU, l-iU > 




8 


9.1232E-3 




9.1489E-3 






2.8701E-5 




10 


9.1174E-3 




9.1621E-3 






4.7919E-5 




12 


9.1029E-3 




9.1587E-3 






5.8863E-5 




14 


9.0883E-3 




9.1502E-3 






6.4679E-5 




16 


9.0776E-3 




9.1404E-3 






6.7453E-5 




extrapolation 


9.0094(4)E-3 




8.9645(5)E-3 




6.64(2)E-5 




conjectured 


*** 




8.9286E-3 




*** 




Table 18 Numerical values of < W2,i-2U > computed 


using the RSOS representation." 


2L 


X = 3 


4 


5 


6 


7 


8 


0.48725 


0.58002 


0.62218 


0.64544 


0.65976 


10 


0.49182 


0.58415 


0.62587 


0.64881 


0.66290 


12 


0.49431 


0.58642 


0.62790 


0.65066 


0.66462 


14 


0.49581 


0.58779 


0.62913 


0.65177 


0.66564 


16 


0.49679 


0.58868 


0.62992 


0.65248 


0.66630 


extrapolation 


0.5001(3) 


0.5917(3) 


0.6326(3) 


0.6547(2) 


0.6683(2) 


conjectured 


0.50000 


0.59161 


0.63246 


0.65465 


0.66815 



In this representation < u, ^2^2 > equals < w, Z_2W2 > and the conjectured value 



of < w2, 1—2^ ^ is equal to W C 



Table 19a Numerical values of < W2, /-2U > computed in the RSOS representation 
with free boundary condition. ° 



2L 


X = 2) 


4 


o 
o 


n A QQ1 o 

u.4oy 


n c:i 1 Qi 


10 


0.45340 


0.52894 


12 


0.46289 


0.54052 


14 


0.46950 


0.54872 


16 


0.47429 


0.55476 


extrapolation 


0.499(3) 


0.582(2) 


conjectured 


0.50000 


0.59161 



° The conjecture value is y^c/2. 



Table 19b Numerical values of < W2, /~2U > computed in the RSOS 
representation with free boundary condition. 



2L 


X = b 


6 


7 


6 


0.51599 


0.53113 


0.54009 


8 


0.54219 


0.55789 


0.56714 


10 


0.56036 


0.57661 


0.58617 


12 


0.57283 


0.58953 


0.59933 


14 


0.59669 


0.61418 


0.62444 


extrapolation 


0.60145(1) 


0.62022(7) 


0.6372(1) 


conjectured 


0.63246 


0.65465 


0.66815 
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